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Abstract 

Algorithms testing two graphs for isomorphism known as yet in com- 
puter science have exponential worst case complexity. In this paper we 
propose an algorithm that has polynomial complexity and constructively 
supplies the evidence that the graph isomorphism lies in P. 

1 Introduction 

Graph isomorphism is a crucial problem in computer science and has been in- 
vestigated intensely in the past. Applications can be found both inside and 
outside the computer science area. Symmetry breaking is an important exam- 
ple, which is a current topic in constraint solving and may be used for instance 
in register allocation, which is a critical performance issue. In chemistry, graph 
isomorphism algorithms are employed for molecule analysis. 
Known algorithms like the ones presented by Nauty [1] and VF [2] do indeed 
indicate an instant isomorphism efficiently in the majority of cases. Yet, for 
each of these techniques, graphs leading to an exponential growth in calculation 
time can be constructed. For a long time, the existence of a method with poly- 
nomial runtime was in doubt, albeit Kobler [3], [4] had expressed that there were 
reasons to believe graph isomorphism not to be NP hard, unlike the subgraph 
isomorphism problem or the travelling salesman problem for which it is assumed 
that fast algorithms do not exist. He pointed out that the problem of finding 
an isomorphism is equally hard to the one of finding all isomorphisms which is 
unusual in NP-hard problems. 

It has been well known [5] that a necessary (but not sufficient) criterion for 
graph isomorphism are identical eigenvalues of the adjacency matrix. This has 
been exploited [6] to develop an algorithm for the class of graphs with bounded 
eigenvalues. 

The matrices will still have same eigenvalues, if you perturbs the diagonal ele- 
ments of the matrices, such that a vertex in the first and in the second graph 
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have same values for the diagonal elements, if they are candidates for an iso- 
morphism. By iterative perturbation of the diagonal elements the two graphs 
can be checked for isomorphism in polynomial time. 

A program for graph isomorphism with polynomial run time was already im- 
plemented by Trofimov and Smolenskii two years ago [7] by using this trick. 

2 Terms and linear algebra 

A graph G = (V, E) consists of a set of vertices V and a set of edges E. Two 

vertices vl and v2 are connected or adjacent iff {vl,v2} G E. 

Two graphs G = (V, E) and G' = (V, E') are isomorph, iff a bijective mapping 

7T : V -» V exists with {vl, i)2}e£o {n(vl), n(v2)} G E' . 

In the adjacency matrix an entry is 1 in the case when the two respective 

vertices are connected through an edge in E, otherwise 0. As we will only 

consider undirectional graphs, the adjacency matrix is symmetric. Kobler [4] 

has proven this to be sufficient as the graph isomorphism problem for directed 

graphs can be reduced to the one for undirected ones. 

A symmetric matrix can be decomposed by means of eigenvalue decomposition 
into A = V * D * V T , D being a diagonal matrix, whose main diagonal entries 
are filled by the eigenvalues of the matrix A. V is a unitary matrix, that is 
V T * V = I. The columns of V are the eigenvectors of A. If V = (t>i, v n ) and 
D = diag(ni,...,n n ), then Aui = faiii. V\,...,v n are orthonormal[8], p. 291, 
thus 



Let {/Ui, ...,/x m } be the set of distinct eigenvalues. Associated with the eigen- 
values in is the eigenspace Si containing the eigenvectors associated with ^ : 
Si = {x G R n | Ax = fiix} . By virtue of the symmetry of A, we have: 

1. If m is an eigenvalue with multiplicity rrn then Si has dimension mj. 

2. The direct sum Si © #2 ... © S m equals R n 

3. If i ^ j then Si and Sj are mutually orthogonal. 

3 Derivation 

If two graphs G and G' are isomorphic, so for the respective adjacency matrices 
A' = P*A*P T must hold, where P denotes a permutation matrix. This entails 
A' k = P*A k * P T for arbitrary k. A' and A have the same eigenvalues, because 
P is an unitary matrix. 

Proposition 1. If A is an adjacency matrix of a graph and P is a permutation 
matrix, then P describes an automorphism of the graph, iff PA = AP. 
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The proof is written in [5]. 

Lemma 1. Let A and A' = P * A * P T be adjacency matrices of isomorphic 
graphs and D a diagonal matrix, so their characteristic polynomials are equal 
and \{A + D) = \{A' + P*D*P T ) holds, i. e. the two matrices have the same 
eigenvalues. 

Proof. A' + P * D * P T = P * (A + D) * P T , applies, which results in the 
eigenvalues and consequently the characteristic polynomials being equal. □ 




Lemma 2. Let A and B be two positive definite matrices, and it furthermore, 
let tt be a permutation tt with Y^j=i(Aji) 2 = Y^j^i-^jm) 2 f or ^ = l,-~ t n> 
then they possess identical eigenvalues and consequently identical characteristic 
polynomials. 

Proof. Let e\, ...,e n be the unit vectors, let fi\, ...,/J, m be the possibly multiple 
eigenvalues and vn, .., Vi mi the eigenvectors belonging to m. 
d can be written as a linear combination of the eigenvectors. 
Let now be e l = Ya=i x ii P v ip> so applies: 

A k * ei = YT=i (/•»? * E™ii xu p vi P ) is the ith column of A k , called af\ We 
observe that 

(fe) T (fe) 

2=1 \ 

I m mi 

=1 p=i 

m m mi n%i' 

= E E * E E Xa * Xil ' P > ' W, V > 
1=1 i'=l p=lp'=l 

m mi m 

= E^ 2 *E*l = E^ 
;=i p=i «=i 

£^ 5 is thereby the component of the vector e, that is within the vector space 
that jii spans. 

Once En are known, we can unambiguously determine the /if from the linear 

system of equations Y^iLiiltfY En = a[ k ^ with i e {1, n}, k G {1, m} 
according to Vandermonde and, as all eigenvalues are positive, the /U; as well. 
Yet, the En still have to be identified. The fact that the eigenvectors form an 
orthonormal basis, yields the condition Y^iLi = 1 for i = 1, to. This prop- 
erty results in a linear system of equations which is unambiguous given that all 
values are positive. 

As this argumentation equally holds for B, both matrices have the same eigen- 
values. □ 
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Gerschgorin's Disc Theorem 1. The eigenvalues of a matrix A — (flij) are 
inside the Gerschgorin discs 

n 

Gi = {z e C : \z - an < ^ Kill »=!,••,« • 

Every disc holds an eigenvalue. 

More information about the theorem can be read, for example, in [9]. 
If A is symmetric, the eigenvalues are real, thus every eigenvalue is in one interval 

[an - E"=i,j^ \ a ij\ > a« + E"=i,^i 

Lemma 3. Let B = D + A with D = diag(d\, d n ) is a diagonal matrix and 
A is an adjacency matrix of an undirected graph. 

Let \di — dj\ > An and di > n — 1 when di ^ dj, otherwise \\B * Ci\\ 2 — 
\\B k * Cj- 1 1 2 for all k. Lf di — dj then j is in the orbit of i, i.e. there is an 
automorphism tjj with j = ip(i) 

Proof. Let N(w) :— {v\v is adjacent to vertex w}, then Bei = di * ej + 
Ei'G-/v(i) e *' • If c^i — dj, the vertices i and j have the same degrees, because 
\\Be l \\ 2 = d\ + £i'eJV(i) I 2 = 4 + Ej'ejvy) l2 = H Se jll 2 - e * can be written 
e i = Y^ij=i x ij v ji where Vi,...,v n are eigenvectors, which are an orthonormal 
basis. Thus 



1 = (e;,e;) = (^XijVj^XijVj) = ^y^XijX U {Vj,Vi) 

3=1 3 = 1 3 = 1 1=1 

n n n 

= X x l ( v i ' v ^ + X x ^ xa ( v i ' v ti = X x % * 1 + X Xi i * = X x2 j 

i=i j=i j=i 



4 



and 



and it is also 



x ij v j > ^ij w j 

\j=l J'=l / 

/ n n \ 

/ n n \ 

\j=l 3=1 I 



n n 



y^y^ nixjjxu (vj,vi) 

3=1 1 = 1 
n 

J2^3 X l = d * 

3=1 



{Be i ) 1 {Be l ) = (Be^Bet) 



\ 3=1 3=1 I 

j n n \ 

tfXijVj ) 

\j=l 3 = 1 I 

In n \ 

y^Xij^ijVj , y^^jXijVj 

\j=l 3 = 1 I 



n n 



= y2yi^3^l X i3 X il l \ V 3^ V l) 
3=1 1=1 

n 

3=1 



where Si is the degree of vertex i. 
Because of Gerschgorin's disc theorem B is positive definite, thus ej and e d - are 
equally distributed to the eigcnspaces as shown in Lemma 2, if di — dj. 
Let E ir = J2je{j\ v es r } x % ^ e * ne ener §y °f e i m the eigenspace S r = {v\Bv = 
/z r u}, then Y^iT=i = 1- We can 1°°^ &t En, Ei m as if it is a probability 
distribution. 

In this case E\p, r ] = Y^T=i ^rEi r =< £i,Bei >= di and Var\p, r ] — E\p%] — 
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EW? = Er=i l4E ir - i VrE ir ] 2 =< BeuBa > - < e u Be, >= df + 5, - 
df = Si. By Chebyshev's inequality there will be Pr(\^i r — di\ > 2n) < < 

= . If an eigenvalue [i r is outside the Gerschgorin disc, then E ir < . 
Let ej and two vectors with "irE ir = Ej ri so the energy of Bei and Bej is also 
equaly distributed to the eigenspaces and so Vr Ei'ew(i) Ei'r = Ej'ew(j) ^j'r- 
For any d let fcj the number of vertices i' G JV(i) with rfj/ = d and kj the number 
of vertices j' G iV(j) with dj/ = d . fcj = kj, because : 
Assumption: kj < hi. 

Let G(d) be the Gerschgorin disk [d — In, + 2n]. 

Let Ec(d)i = ^2re{r\ fJ , r eG(D)} Eri the energy in the eigenspaces with eigenvalues 
inside the Gerschgorin disc, then 

E G(d^ > - ^) > *i - \ 

i'eJV(i) 

and 

X! <kj + {N{j)-kj)^ <kj + \ . 

j'eN(j) 

From that follows fcj = fcj. 

So there is a bijection between adjacent vertices of i and j with Z and -0(0 with 
di = rf^(i) • This way we find a bijection on every vertex of the graph, which is 
described by a permutation matrix P. If i and j are vertices with ej = Pei, 
then Bu — di = dj = Bjj and so Ei r = Ej r for every r G {l,...,m}. Let 
e, = E^Li an< i e i — E^Li u j>, where Ui r and Uj r are eigenvectors in the 
eigenspace S r , then ||ui r || 2 = E ir = Ej r = ||u ir || 2 . 

Vr Uj r = Pui r , because A k P^2 r Hr u jr = A k Pei = A k ej = E r l^r u jr f° r every 
k and \i r > are distinct eigenvalues, thus a linear system of equations with 
exactly one solution is build. 

This Causes PBei — f*ErLl MrUir = E^Ll Mr-P* M ir = E^Ll Hr-Ujr = B*ej = 

B * Pej, thus PB = BP. For the adjacency matrix A = B — diag{d\, d n ) we 
get PA = AP and so P is an automorphism, as proven in [5] . 

□ 

4 Algorithms 
4.1 Algorithm 1 

input : adjacency matrices of graphs A und A' 
output: graphs isomorphic ? 

<71i := {l,..,n} 
C2i := {l,..,n} 
Al := A 
A2 := A' 
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repeat 

c := number of C'l x 
for i:=l to n 

Al[i,i] := 4n * x, where i G Cl x 
A2[i,i] := 4n * x, where i G C2 X 
for fc:=l to n 
for i:=l to n 

al^:= ith column vector of Al k 
a2^:= ith column vector of A2 k 

Al\m := alf )T alf 

A2[i\[k] := a2f )T a2f ) 
sort lists ^41 and A2 lexicographically 
compare Al and A2 
if sorted lists are not equal 
then return false 

put vertices into equivalence classes Cl x , such that i = j iff Al[i] = Al[j] 
put vertices into equivalence classes C2 X , such that i = j iff A2[i] = A2[j] 
and Al[l] = A2[i] iff I G Cl^ and i G C2,; 

until c = number of C x 

return true 

If the graphs are isomorph, the algorithm will return true, because not only 
the adjacent matrices have the same spectra, also the matrices Al ad A2 have 
the same spectra in every iteration step because of Lemma 1. 
If Algorithm 1 returns true, we can look at the graph as two not connected 
parts of a big graph. Because of Lemma 3 we find for every vertex in the first 
graph a vertex in the second graph, with is in the same orbit, belonging to the 
big graph. That way, the algorithm construct an isomorphism. 

4.2 Algorithm 2 

Algorithm to construct an isomorphism 

input : adjacency matrices of two isomorphic graphs A und A' 
output: graph isomorphism, if exists 

make A1,A2,C1,C2 and c as described in Algorithm 1 
while c < n 

get x with |C1 X | = \C2 X \ > 1 

get an i G Cl x and j G C2 X 

Al[i,i] := A2\j,j] := 4 * (c + 1) * n 

for k:—l to n 
for i:=l to n 
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al\ := ith column vector of Al k 
a2^:= ith column vector of A2 k 

Ai\m ■■= «if )T a if } 

A2[i][k] := a2f )T a2f ] 

put vertices into equivalence classes Cl x , such that i = j iff = .41 [j] 
put vertices into equivalence classes C2 X , such that i = j iff .42 [i] = .42 [j] 

and il[Z] = i2[i] iff Z e Cl x and i g C2 X 
c := number of C'l x 

sort lists Al and A2 lexicographically 

return permutation induced by the sorted lists 

When c = n, we have only single eigenvalues for Al and A2. In this case, 
there is only one solution for Al = P*A2*P T with P unitary, the permutation 
the algorithm has produced. Now A = P * A' * P T because of Lemma 1 and so 
the permutation is an isomorphism of the graphs. 

5 Complexity 

The main effort of Algorithm 1 consists in n matrix multiplications in every 
iteration of the outer loop. The outer loop is a fixed point iteration, because 
it terminates, when no more equivalence classes can be constructed. In case of 
constructing new classes, the number of classes will grow. Algorithm 1 must 
not make more then n iterations in the outer loop. 

In the outer loop of Algorithm 2 vertex i of the first graph and vertex j of 
the second graph will put into a single class separated to the other vertices. 
Algorithm 2 also have not more then n iterations in outer loop. 
Both algorithms consists in n 2 matrix multiplications in n iterations of the outer 
loop. Thus we have to perform 0(n 8 ) integer additions or multiplications. 
The maximum eigenvalue of Al is according to Gerschgorin's circle theorem 
< 4 "(" 2 t ~ 1 ) = 2n 2 + 2n. Consequently the maximum eigenvalue of Al n < (2n 2 + 
2n) n . Since all matrix entries are positive, none of them can be bigger than the 
maximum eigenvalue and therefore < (2n 2 + 2n)™. The maximum value occurs 
in calculating the sum of the quadratic entries and is < n * (2n 2 + 2n) 2n . The 
coding length for an entry is therefore (6n) * (log(n) + C) = 0(n * log(n)). 
Algorithm 2 needs also at most 0(n 2 ) matrix multiplications with maximum 
coding length of an entry in 0(n* log(n)). 

6 Literature 

1. McKay, Brendan D. (1981), Practical graph isomorphism, Congressus Nu- 
merantium 30: 10th. Manitoba Conference on Numerical Mathematics 
and Computing (Winnipeg, 1980) 



8 



2. amalfi.dis.unina.it/graph/db/papers/benchmark.pdf 

3. Kobler, Johannes; Schoning, Uwe & Toran, Jacobo (1992), Graph isomor- 
phism is low for PP, Computational Complexity 2 (4) 

4. Kobler, Johannes; Schoning, Uwe & Toran, Jacobo (1993), The Graph 
Isomorphism Problem: Its Structural Complexity, Birkhauser, ISBN 978- 
0817636807. 

5. Biggs, Norman (1974), Algebraic Graph Theory, Cambridge University 
Press 

6. Babai, Laszlo; Grigoryev, D. Yu. & Mount, David M. (1982), Isomorphism 
of graphs with bounded eigenvalue multiplicity, pp. 310-324, Proceedings 
of the 14th Annual ACM Symposium on Theory of Computing 

7. Trofimov, M. I.; Smolenskii, E. A. (2005), Application of the electronega- 
tivity indices of organic molecules to tasks of chemical informatics, Russian 
Chemical Bulletin, International Edition, Vol. 54, No. 9, pp. 2235-2246 

8. Bronstein, I. N.; Scmendjajew, K. A.; Musiol, G.; Muhlig, H. (2001), 
Taschenbuch der Mathcmatik, 5. Auflage, Verlag Harri Deutsch, Thun 
und Frankfurth am Main 

9. Plato, Robert (2003), Concise Numerical Mathematics, Graduate Studies 
in Mathematics, Volume 57, American Mathematical Society 



9 



